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A new method for a precision measurement of the CKM matrix element \V u b\ is discussed, which 
combines good theoretical control with high efficiency and a powerful discrimination against charm 
background. The resulting combined theoretical uncertainty on \V u b\ is estimated to be 10%. 
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I. INTRODUCTION 



One of the main goals of B physics is to perform a 
detailed study of the flavor sector of the electroweak the- 
ory This can be achieved by accurate measurements of 
the parameters of the unitarity triangle. The angle f3 
has already been determined precisely at the B factories. 
Because \V u b\ is the side opposite to this angle, an ac- 
curate measurement of this parameter will directly test 
the Standard-Model picture of CP violation. It is there- 
fore desirable to device methods for extracting \ V u b\ with 
theoretical uncertainties at the 5-10% level. 

Inclusive semileptonic B — > X u I _ v decays offer the 
cleanest determination of | V u b \ ■ In order to discriminate 
against the large charm background, one has to apply 
tight experimental cuts. So far, cuts on the charged- 
lepton energy Ei > (M§ — M]^)/{2Mb), the hadronic 
invariant-mass squared sh < Mp,, and the dilepton mass 
squared q 2 > {Mb — Mr,) 2 have been employed. Only the 
Ei cut can be applied without neutrino reconstruction. 
Unfortunately, it has a very low efficiency and is there- 
fore theoretically disfavored. The hadronic-mass cut is in 
principle the ideal separator between B — > X u l~v and 
B — > X c l~v events. However, in practice one has to 
lower the cut due to the experimental resolution on sh, 
thereby reducing the efficiency. In this Letter we elab- 
orate on a highly efficient method for a precision mea- 
surement of \V u b\ recently proposed in [1], which offers 
several advantages over the existing strategies. 

In the relevant region of phase space, differential B — > 
X u l~v decay rates can be factorized into perturbatively 
calculable functions H L and Ji , and universal shape func- 
tions Si [2-4]. The hard functions Hi account for quan- 
tum fluctuations at high scales fih ~ mb, whereas the 
jet functions Ji describe the properties of the final-state 
hadronic jet at scales Hi 
Hi 



'sh ~ ^m b AQ C D- Both 
and Ji are known to one-loop order in perturba- 
tion theory [1,5]. The functions Si parametrize the non- 
perturbative physics below the intermediate scale \Xi. 
In [1] we have derived several model-independent con- 
straints on the form of the leading-order shape function 
S. This significantly reduces the hadronic uncertainties 
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FIG. 1. Phase space for the variables P_ and P+. The dif- 
ferent shadings separate the regions where sh = P+P- < M% 
(dark gray) and sh > Md (light gray). The two solid lines 
correspond to the exclusive decays B — » D^*'l~D, whereas 
the dotted line displays the lower boundary for inclusive 
B — > X c l~v decays as predicted by the heavy-quark expan- 
sion. The dashed line corresponds to P+ = M^/Mb- 

in predictions for decay spectra and event fractions. In 
principle, the shape function can be constrained further 
using information about the B — > X s "f photon spectrum. 

The hadronic physics in B — > X u l~ v is most naturally 
described using the variables P± = Eh T \Ph\, where 
Eh and Pjj are energy and three-momentum of the final- 
state hadronic jet. The phase space in these variables is 
depicted in Figure 1 (for P + < 3 GeV). The vast major- 
ity of events is located in the "shape-function region" of 
small P + and large P_, corresponding to Eh 3> ^Jsh- 
In this region a rigorous theoretical description can be 
obtained using soft-collinear effective theory [6], which 
provides a systematic expansion in the small parameter 
Aqcd/w;,. All B — > X c l~v events are confined to the 
light-shaded region. Cutting on the variable P + allows 
one to isolate the shape-function region from the rest of 
phase space. Choosing P+ < Mjj/Mb ~ 0.66 GeV elimi- 
nates the charm background entirely while maintaining a 
high efficiency. This cut samples the same hadronic phase 
space as the lepton-endpoint cut Ei > 2.31 GeV, but it 
contains significantly more events. Here we propose to 
employ such a cut to determine \V u b\- 
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The hadronic-mass cut sh < M|, also fully contains 
the shape-function region. However, in addition it in- 
cludes a triangular region of larger P + . Since the contour 
Sh = borders the charm background, the observed 
event fraction is very sensitive to "leakage" of charmed 
events due to imperfect experimental resolution. The P+ 
cut, on the other hand, provides a "buffer zone" against 
charm background. 



II. CALCULATION OF THE EVENT FRACTION 

The fraction Fp of all B — > X u l~v events with 
hadronic light-cone momentum P+ < Ap is given by 



Fp(Ap) = T(a) e VH ^ h '^ 



m b\ a [ Ap 
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and the functions T, Vp, Hi, and Jj can be found 
in [1]. This result is valid in the shape- function region 
Ap ~ Aqcd and at leading power in the heavy-quark ex- 
pansion. Sublcading contributions can be calculated sys- 
tematically using the techniques developed in [1]. Large 
logarithms have been resummed into a leading-order Su- 
dakov factor T(a) e Vn and next-to-leading order func- 
tions Hi and Jj. The result for Fp is formally indepen- 
dent of the two matching scales ph and p^. 

The fraction Fm of events with hadronic invariant 
mass such that sp < so receives a contribution from the 
box-shaped region below the dashed line in Figure 1 as 
well as the triangular dark-shaded region above that line. 
The former simply yields Fp(A s ) with A s = so/Mp, 
while the latter involves a weighted integral over the 
shape function up to ^/so [1]. Unlike for the P + spec- 
trum, for which knowledge of the shape function in the 
range < ui < Ap <~ Aqcd is needed, here the functional 
form of the shape function is required over a wider range. 
Because the collinear expansion breaks down near the tip 
P + ~ P_ ~ yfs~o of the triangular region, it is an open 
challenge how to construct a systematic heavy-quark ex- 
pansion for Fm(so)- 

Our prediction for the fraction of events with the op- 
timal cut P + < M 2 D jM B is 



F P = (79.6 ± 10.8 ± 6.2 ± 8.0)% , 



(3) 



where the errors represent the sensitivity to the shape 
function, an estimate of 0(a 2 ) contributions, and power 
corrections, respectively. The uncertainty due to our ig- 
norance of the shape function is obtained by varying the 



parameters A = Mp — m b and /i^. determining its first 
two moments (5F P — lg' 2 %), as well as the functional 
form of S (SFp = Here, m b and are de- 

fined in the shape- function scheme [1]. Higher-order per- 
turbative effects are estimated by studying the depen- 
dence on the matching scales Hi and jit/j, varied in the 
ranges 1.25 GeV < m < 1.75 GeV (SF P = t£g%) and 
m fc /V2 < [i h < y/2m b (8F P = The pertur- 

bative uncertainty could be reduced by computing the 
0(a 2 s ) corrections to Eq. (1). Power corrections will be 
discussed in more detail below. 

The CKM-matrix element \V u b\ can be extracted by 
comparing a measurement of the partial rate T U (P + < 
Ap) with a theoretical prediction for the product of the 
event fraction Fp and the total inclusive B — > X u l~v 
rate. The resulting theoretical uncertainty on \V u b\ is 



S\V U , 

\v ub 



= (±7 ± 4 ± 5 ± 4)% , 



(4) 



where the last error comes from the uncertainty in the 
total rate [7,8] . Because of the large efficiency of the P + 
cut, weak annihilation effects [9] have an influence on 
\V u b\ of less than 2% and can be safely neglected. 

At leading power in Aqcd/w6 the shape- function un- 
certainty could be eliminated by relating the P+ spec- 
trum to the B — > X s "/ photon spectrum, both of which 
are given (at tree level) directly by the shape function. 
A prototype of such a relation can be found in [1] , where 
the charged-lepton energy spectrum was expressed as a 
weighted integral over the P + spectrum. Using similar 
methods, it will be possible to construct a shape-function 
independent relation of the form 
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(5) 



where at tree level the weight function is simply 
w(A,Ej) = 1. At next-to-leading order the energy de- 
pendence of this function resides in a single logarithm 
ct s (ni) ln[mf,(A — Mp + 2E~ f )/^i 2 ]. A similar expression 
relating the hadronic-mass spectrum with the B — > X s -f 
photon spectrum would be far more complicated and re- 
quire input of the photon spectrum beyond the region 
where it is currently experimentally accessible. 

In order to study power corrections it is instructive 
to resort to a local operator-product expansion (OPE), 
which is valid for large values of P + and P_ . The result- 
ing prediction for the normalized B — > X u l~v spectrum 
in the variable p + = p + /mb = (P+ — A)/m& reads 
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We include the contributions from dimension-3 operators 
at 0(a s ) and those from dimcnsion-5 operators (whose 
matrix elements are proportional to the heavy-quark ef- 
fective theory parameters Ai^) at tree level, using [10,11]. 
If the p + spectrum is integrated without a weight func- 
tion, the tree-level power corrections from dimension-5 
operators vanish. This is in accordance with the fact that 
subleading shape functions have zero norm at tree level 
[12]. For Aqcd < Ap « m;, we obtain with A = A P -A 
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For Ap ~ 1 GeV, the power correction in the second line 
leads to an enhancement of the fraction Fp by 5-10%. 
This is of similar magnitude as tree-level estimates of 
(zero-norm) subleading shape-function effects on the E\ 
and sh spectra [13,14]. 

The leading-order term in Eq. (8) can also be obtained 
from Eq. (1) by taking the limit Ap 3> Aqcd- Interest- 
ingly, such an analysis uncovers that in this kinematic 
range there is an enhanced class of power corrections of 
the form (Aqcd/ A)" with n > 2, which arise first at 
0(a s ) (see also [5]). The leading corrections to the ex- 
pression above are given by 



5F P (Ap) = -^-[ 2 Ing- 



raft 



A. 

3A 2 



(9) 



It follows that for sufficiently large A the event fraction 
can be expressed as a double expansion in A/ m b and 
Aqcd/ A. Such a multi-scale OPE is non-trivial and left 
for future work. 



III. CHARM BACKGROUND 

One of the main advantages of the P+ spectrum over 
the hadronic-mass spectrum is a better control of the 
charm background. In order to study this in more detail, 
we investigate the OPE prediction for the normalized p + 
spectrum in B — > X c l~D decays, 
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dT c = 2{g-p 2 + , 

df>+ f{g)p+ 

+ gp+(3 - 8p+ + 3p 2 + ) - g 2 (2 - 3p+) 



(10) 



where g < p+ < yfg with g — (m c /m b ) 2 , and 

f(g) = l-8g + 8g 3 - g 4 -12g 2 lng. (11) 

We only include tree-level contributions from dimension- 
3 operators. In this approximation, inclusive charm 
events are located along the dotted line in Figure 1. Be- 
cause 0(a s ) corrections redistribute these events into the 
light-gray segment above that line, the integral over the 
spectrum in Eq. (10) serves as an upper bound on the in- 
clusive charm background. In the same approximation, 
the hadronic mass distribution is given by 
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Here, s = s#/to 2 , e = A /nib, and phase space is such 
that (y/g + e) 2 < s < (g + e)(l + e). 

The results for the event fractions Fp and Fm are 
summarized in Figure 2. The bands show the inclusive 
B — > X u l~v event fractions for different shape- function 
models [1], while the solid lines give predictions for the 
inclusive B — > X c l~v background spectra as obtained 
from Eqs. (10) and (12). The latter are normalized 
to their total decay rate T c , which is about 60 times 
larger than r„. In the case of the fraction Fm, the 
inclusive charm background starts right at the thresh- 
old A s = (m c + A) 2 1 M B =± M^/Mp. In the case of 
the Fp event fraction, on the contrary, the inclusive 
charm background starts at a value Ap = m 2 c /m b + A w 
0.96 GeV, which is significantly larger than the value 
Ap = M|, I Mp ~ 0.66 GeV above which final states con- 
taining charm mesons arc kincmatically allowed. This is 
a reflection of the fact that there is a gap between the 
dashed line representing the P + cut and the dotted curve 
representing the onset of the continuum of charmed final 
states in Figure 1. In reality this gap is filled by exclu- 
sive modes containing the hadronic final states D, D* , or 
Dtt, Dtttt. The exclusive contributions to the hadronic 
P + spectrum from B — > D^*H~9 decays are 
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where M 2 (t) /M B < P+ < M DM and r {sf) = M DM /M B - 
The recoil variable w — v ■ v' is given by M^ )(t) + P 2 = 
2wM D (,)P+. The corresponding contributions to the 
hadronic mass spectrum are given by 



ds H 



= T DM ■ 5{s H - M D( ,)) . 



(14) 
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FIG. 2. Fraction of B — > X u l v events with P+ < Ap (left) and sh < so (right). The bands correspond to the values 

A = 0.63 GeV (dark), 0.70 GeV (medium), and 0.56 GeV (light). Their width reflects the sensitivity to the value of /z 2 . 

varied between 0.20 and 0.34 GeV 2 . The curves represent the background from inclusive B — > X c l~P (solid) and exclusive 

B — * D^l~V decays (dashed), normalized to the total inclusive semileptonic charm rate. 



The exclusive contributions to the event fractions are 
given by the dashed lines in Figure 2. To obtain these 
curves we have used an ansatz for the form factors 
T D (*) (w) that is consistent with experimental data on 
the recoil spectra and branching fractions. The fact that 
these exclusive modes are very well understood should 
help to model the background. The smooth onset of the 
Z)(*) background is a direct consequence of the fact that 
the ideal P+ cut touches the charm region at only a sin- 
gle point in phase space (see Figure 1). On the contrary, 
the region of phase space when applying the sh cut bor- 
ders the background along the curve separating the light - 
and dark-shaded regions, which leads to a step increment 
in the event fraction Fm- As a consequence, one needs 
to move away from the ideal cut sh — M|, because of 
smearing effects due to experimental resolution. It is our 
hope that it will be possible to stay closer to the ideal cut 
when performing a P+ analysis. In this case, both the 
P + and sh discrimination methods lead to comparable 
efficiencies and shape-function uncertainties. 

IV. CONCLUSIONS 

We have presented a new method for a precision mea- 
surement of the CKM matrix element \V u b\ in inclusive 
B — > X u l~v decays, using a cut on the hadronic light- 
cone variable P+ = Eh — \Ph\- This method has been 
compared to a strategy using a cut on hadronic invari- 
ant mass. While both approaches are highly efficient 
and comparable in their sensitivity to non-perturbative 
physics, a measurement of the P + spectrum is favored 
in view of its theoretical potential. The reasons arc 
twofold: First, it offers a simpler construction of shape- 
function independent relations to the photon spectrum 
in B — > X s -y, thus allowing for a theoretically clean de- 
termination of \V u b\. Secondly, the P+ spectrum can be 
evaluated within a systematic framework, which makes 
the calculation of power corrections to Eq. (1) feasible. 



This will be necessary in order to lower the theoretical 
uncertainty on \V u b\ to the 5% level. Such an ambitious 
goal might also require to include at least part of the 
0(a 2 s ) corrections. 
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